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I. INTRODUCTION 



This conference is devoted to the physics of waves moving through a medium which 
affects them as would a background metric. There is therefore a natural analogy with waves 
propagating in a gravitational field — but we should take care not to push the analogy 
too far. These systems differ from real general relativity in that the medium constitutes 
a physical aether. While we may ignore the aether for many purposes, occasionally it is 
important. For example, if we wish to compute forces exerted by the waves, we must take 
into account any stress transmitted by the background medium. 

The natural tool for computing forces by tracking the flux of energy and momentum 
is the energy-momentum tensor. This is best defined as the functional derivative of the 
action with respect to the background metric — but we have two metrics at our disposal: 
the spacetime metric and the acoustic or other metric which we are exploiting for our GR 
analogy. Despite the temptation to believe otherwise, we should remember that it is only 
by differentiating with respect to the "real" metric that we obtain "real" energy and "real" 
momentum. When we differentiate with respect to the analogy metric, we obtain the density 
and flux of other quantities. These are usually the pseudoenergy and the pseudomomentum 

E- 

Failure to distinguish between real energy and momentum and pseudoenergy and pseu- 
domomentum has caused much confusion and controversy over the years. Consider the 
formula for the momentum of a photon in a dielectric: should the refractive index go in the 
numerator, where it was placed by Minkowski, or in the denominator, as argued by Abra- 
hams? This dispute was not resolved until Blount identified Minkowski's expression with 
the pseudomomentum and Abrahams' with the true momentum, including the mechanical 
momentum of the dielectric ||3|||] . 

Although I will not address the problem here, my initial motivation for thinking about 
these topics was the fear that a similar confusion lies behind some recent controversies 
involving the Iordanskii force. This force, which has an appealing GR analogue in the 
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gravitational Aharonov-Bohm effect of a spinning cosmic string ||, is supposed to act on a 
vortex in a superfluid when it moves relative to the normal component of the fluid. 

A related issue, and one that I will address, lies at the heart of the two-fluid model for 
a superfluid or Bose condensate. It is one of the fundamental assumptions of the two-fluid 
model that the phonons in a fluid possess momentum Kk, and that, unlike that of phonons 
in a solid, this momentum is true Newtonian momentum, mv. This assumption is essential 
because we wish to identify the phonon momentum density with the mass current, both 
being equal to pv. The desired identification is supported by the approximate solution of 
Bogoliubov's weakly interacting Bose gas model, in which the phonon creation operator, a^, 
appears to be create bona-fide momentum, so it is quite unnerving to discover that in the 
literature of fluid mechanics the attribution of real momentum to a sound wave is regarded 
as a naive and dangerous fallacy. A particularly forceful statement of this opinion is to be 
found in the paper [|7| On the Wave Momentum Myth, by Michael Mclntyre. 

In the present article we will focus the difference between the true momentum and the 
pseudomomentum associated with the acoustic metric. The first part is a general pedagogical 
account of the distinction between momentum and pseudomomentum, and the circumstances 
under which the latter may be used for computing forces. The second part will discuss the 
energy and momentum associated with sound waves in a background flow ||. 



The distinction between true momentum and pseudomomentum is especially clear when 
we consider the problem of transverse vibrations on an elastic string. The action 



II. MOMENTUM AND PSEUDOMOMENTUM 




(2.1) 



gives rise to the familiar wave equation 



py - Ty 
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(2.2) 
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with c = yjT/p being the wave speed. By manipulating the wave equation we can establish 
two local conservation laws. The first, 



|{i^1 + J^™^ < 2 - 3 > 



is immediately recognizable as an energy conservation law, with the flux —Tyy' being the 
rate of doing work by an element of the string on its neighbor. The second, 



is slightly more obscure in its interpretation. 

In a relativistic system the appearance of —pyy' = (energy flux)/c 2 as the local density 
of a conserved quantity would not be surprising. The symmetry of the energy-momentum 
tensor requires that T j0 = T 0j , so the energy flux T j0 is also (after division by c 2 ) the density 
of 3-momentum. Here, however, we are dealing with non-relativistic classical mechanics - 
and with transverse waves. Whatever the quantity —pyy' may be, it is not the density of 
the x component of the string's momentum. It is instead the density of pseudomomentum. 

To understand the origin of pseudomomentum, observe that our elastic string may be 
subjected to two quite distinct operations either of which might be called "translation in 
the x direction" : 

• An operation where the string, together with any disturbance on it, is translated in 
the x direction. 

• An operation where the string itself is left fixed, but the disturbance is translated in 
the x direction. 

The first operation leaves the action invariant provided space is homogeneous. The 
associated conserved quantity is true Newtonian momentum. The second operation is a 
symmetry only when both the background space and the string are homogeneous. The 
conserved quantity here is pseudomomentum. Such a distinction between true and pseudo- 
momentum is necessary whenever a medium (or aether) is involved. 




(2.4) 
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Adding to the confusion is that, although the pseudo momentum is conceptually distinct 
from the the true momentum, there are many circumstances in which changes in pseudo- 
momentum can be used to compute real forces. As an example consider a high speed train 
picking up its electrical power from an overhead line. 

F 
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Fig 1. A high-speed train . 

The locomotive is travelling at speed U and the pantograph pickup is exerting a constant 
vertical force F on the power line. We make the usual small amplitude approximations 
and assume (not unrealistically) that the line is supported in such a way that its vertical 
displacement obeys an inhomogeneous Klein- Gordon equation 



py - Ty" + pQ 2 y = F5(x - Ut), 



(2.5) 



with c = yT/p, the velocity of propagation of short wavelength disturbances. 

If < c, the vertical displacement relaxes symmetrically about the point of contact. 
Once U exceeds c, however, the character of the problem changes from elliptic to hyper- 
bolic, and an oscillatory "wake" forms behind the pantograph. As with all such wakes, the 
disturbance is stationary when viewed from the frame of the train. With this in mind, we 
seek a solution to ( f2.5| ) of the form y = y(x — Ut). Since the overhead line is undisturbed 
ahead of the locomotive, we find 

'"y Fc Q'y 
y = — — — sin — (Ut — x), x < Ut 

y = 0, x>Ut. (2.6) 

Here 7 = (U 2 /c 2 — 1)~^ is the Lorentz contraction factor modified for tachyonic motion. The 
condition that the phase velocity, u/k, of the wave constituting the wake be equal to the 



forward velocity of the object creating it is analogous to the Landau criterion determining 
the critical velocity of a superfluid. There are no waves satisfying this condition when U < c, 
but they exist for all U > c. 

In the wake, the time and space averaged energy density {£) = {\py 2 + \Ty' 2 + \Q 2 y 2 ) 
is given by 

<0 = ^W^V (2.7) 



T 1 \T, 

The expression for the pseudomomentum density for the Klein-Gordon equation is the same 
as that for the wave equation, and the average pseudomomentum density is 

1 4 /F^ 2 



(-pyy') = ^Pl\j;) U. (2.8) 

Because energy is being transfered from the locomotive to the overhead line, it is clear 
that there must be some induced drag, on the locomotive. This is most easily computed 
from energy conservation. The rate of working by the locomotive, F d U, must equal the 
energy density times the rate of change of the length of the wave-train. Thus 

F d U=(E)(U-U g ) (2.9) 

where 

a 5 = f£ = ^ = £ (2.10) 
dk uj U 

is the group velocity of the waves, and so the speed of the trailing end of the wake wave-train. 
After a short calculation we find that the wave-induced drag force is 

Since the average pseudomomentum density turned out to be the average energy density 
divided by U, we immediately verify that we get exactly the same answer for F d if we equate 
the wave drag to the time rate of change of the total pseudomomentum. 

For the sceptic we note that we may also obtain the same answer by a more direct 
evaluation of the force required to deflect the overhead wire. 
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Fig 2. Force on overhead wire. 

From the solution fl2.6|) we see that the force F is related to the angle of upward deflection, 
6, by 9 m tanO = —y'(Ut,t) = ^y 2 F/T. By balancing the acceleration of the power line 
against the force and the tension in the line, we see that the force cannot be exactly vertical, 
but must be symmetrically disposed with respect to the horizontal and deflected parts of the 
line. The force exerted by the pantograph thus has a small horizontal component .F sin 0/2. 
The wave drag is therefore F^ = |F 2 7 2 /T, as found earlier. 

Since there is a real horizontal force acting on the wire, true Newtonian momentum must 
also be being transfered to the wire. Indeed the wire behind the train is being stretched, 
while that in front is being compressed. A section of length 2ci ong (t — 1 ), where ci ong is the 
velocity of longitudinal waves on the wire and (t — t ) is the elapsed time, is in uniform 
motion in the x direction. Since usually ci on g ^> c, this true momentum is accounted for by 
an almost infinitesimal motion over a large region of the wire. The pseudomomentum and 
the true momentum are to be found in quite different places — but the divergence of their 
flux tensors, and hence the associated forces, are equal. 

III. RADIATION PRESSURE 

In this section we will consider the "radiation pressure" exerted by sound waves incident 
on an object immersed in the medium. This is a subject with a long history of controversy 
PJT0|. The confusion began with the great Lord Rayleigh who gave several inequivalent 
answers to the problem. Our discussion will follow that of Leon Brillouin JO]] who greatly 
clarified the matter. We begin by considering some analogous situations where the force is 
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exerted by transverse waves on a string. 

Consider a standing wave on a semi-infinite elastic string. 




Fig 3. A vibrating string exerts a force on a bead. 

We have restricted the vibration to the finite interval [0, L\ by means of a frictionless 
bead which forces the transverse displacement of the string to be zero at x — L, but allows 
free passage to longitudinal motion, and so does not affect the tension. 

Suppose the transverse displacement is 



y = A sin cut sin 



irnx 



(3.1) 



with uj = cKnjL. The total energy of the motion is 

E = jdx^-py* + l -Ty»} = \p^A*L. 



(3.2) 



If we alter the size of the vibrating region by slowly moving the bead, we will alter the 
energy in the oscillations. This change in energy may be found by exploiting the Boltzmann- 
Ehrenfest principle which says that during an adiabatic variation of the parameters of a 
harmonic oscillator the quantity E/uj remains constant. Thus 5(E/uj) = or 5E = (J^j 5u. 
To apply this to the string, we note that 

5L 



5u = 5 



cnn 



—uj- 



L J L 

The change in energy, and hence the work we must do to move the bead, is then 



(3.3) 



5E 



(3.4) 



The "radiation pressure" is therefore E/L = (S), the mean energy density. This calculation 
can be confirmed by examining the forces on the bead along the lines of Fig. 2. 
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The average density of pseudomomentum in each of the two travelling wave components 



bounce off the bead. Thinking through this example shows why pseudomomentum can be 
used to compute real forces: On a homogeneous string the act of translating the bead and 
the wave together, while keeping the string fixed, leaves the action invariant. The associated 
conserved quantity is the sum of the pseudomomentum of the wave and the true momentum 
of the bead. 

Keeping track of pseudomomentum cannot account for all forces, however. There is 
another plausible way of defining the radiation pressure. This time we use a finite string 
and attach its right-hand end to a movable wall 



Now, as we alter the length of the string, we will change its tension, and so alter the 
value of c. We must take the effect of this into account in the variation of the frequency 



of the standing wave is ±|(£)/c. The radiation "pressure" can therefore be accounted for 
by the 2 x c x |(£)/c ra ^ e °f change of the pseudomomentum in the travelling waves as they 




5L 



L 



Fig 4. Another way to define radiation pressure. 




(3.5) 



The change in the energy of the vibrating system is therefore 




(3.6) 



The radiation pressure is thus given by 




(3.7) 



This force is in addition to the steady pull from the static tension T in the string. 
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The generality of the Boltzmann-Ehrenfest principle allows us to apply the previous 
discussion with virtually no changes to compute forces exerted by a sound wave. We need 
no explicit details of the wave motion beyond it being harmonic. The two ways of defining the 
radiation pressure for a vibrating string correspond to two different experimental conditions 
that we might use for measuring the radiation pressure for sound waves. The movable end 
condition corresponds to what is called the Rayleigh sound pressure. 



Fig 4. The Rayleigh Radiation Pressure. 
We establish a standing wave in a cylinder closed at one end and having a movable piston 
at the other. Moving the piston confining the sound wave changes both the wavelength of 
the sound and the speed of propagation, producing a sound radiation pressure 



The subscript S on the derivative indicates that it is being taken at fixed entropy. As with 
the string, this radiation pressure is in addition to the equilibrium hydrostatic pressure, Po, 
on the piston. 

The analogue of the string with the sliding bead leads to the Langevin definition of the 
radiation pressure. Here we insert a bypass so that moving the piston confining the sound 
wave does not change the density or pressure of the fluid. The radiation pressure on the 
piston is simply (£). 





(3.8) 




Fig 5. Langevin Radiation Pressure. 
The difference in the two definitions of radiation pressure arises because, if we keep the 
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mean pressure fixed, the presence of a sound wave produces an 0(A 2 ) change in the volume 
of the fluid. If, instead, the mean density is held fixed, as it is in the Rayleigh definition, 
then this volume change is resisted by an additional hydrostatic pressure on the walls of the 
container. 

The most common way of measuring sound pressure involves a sound beam in an open 
tank of fluid. Since the fluid is free to expand, this corresponds to the Langevin pressure. 



Fig 6. The usual experimental situation. 
The radiation pressure is, in reality, a radiation stress 

The anisotropic part depends on the wave-vector k of the sound beam, and may be ac- 
counted for by keeping track of the pseudomomentum changes. The isotropic part cannot 
be computed from the linearized sound-wave equation since it requires more information 
about the equation of state of the fluid medium than is used in deriving the wave equation. 
The extra information is encapsulated in the parameter (fjjj^;) , a fluid-state analogue of 
the Gruneisen parameter which characterizes the thermal expansion of a solid. When the 
experimental situation is such that this isotropic pressure is important, the force associated 
with the sound field cannot be obtained from the pseudomomentum alone. 

IV. MASS FLOW AND THE STOKES DRIFT 



Further confusion involving momentum and pseudomomentum in acoustics is generated 
by the need to distinguish between the Euler (velocity field at a particular point) description 
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of fluid flow, and the Lagrangian (following the particles) description. 
Suppose the velocity field in a sound wave is 

t) — Acos(kx — ut). (4.1) 

Using the continuity equation d x pv + d t p = 0, setting p = p + pi, and approximating 
pv ~ PoV(i), we find that 

Pl = - Po kA cos(kx - cut) + 0(A 2 ). (4.2) 

CU 

The time average of the momentum density pv is therefore 

(pv) = (Piv {1) ) = k -^\A\ (4.3) 

CU z 

to 0(A 2 ) accuracy. This Newtonian momentum density is clearly non-zero, and numerically 
equal to the pseudomomentum density. Here, unlike the case of the elastic string, there is 
only one velocity of wave propagation and so the pseudomomentum and true momentum, 
although logically distinct, are to be found in the same place. 

Further (pv) is both the momentum density and the mass-current. A nonzero average for 
the former therefore implies a steady drift of particles in the direction of wave propagation, 
in addition to the back-and-forth motion in the wave. We can confirm this by translating 
the Eulerian velocity field v^ = Acos(kx — cut) into Lagrangian language. The trajectory 
£ (t) of a particle intially at xq is the solution of the equation 

f = v {1) (Z(t),t) = Acos(kt - cut), m = x . (4.4) 

Since the quantity £ appears both in the derivative and in the cosine, this is a nonlinear 
equation. We solve it perturbatively by setting 

£(t) = x + A E,(t) + A 2 E 2 (t) H . (4.5) 

We find that 

Si(i) = sin(fcco — cot). (4.6) 

cu 
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Substituting this into (|4.4|) we find 

h 

sin 2 (/cx - cot). (4.7) 



dE 2 k . 2 



dt u 

Thus S2 has a non-vanishing time average, k/(2u), leading to a secular drift velocity vl = 
^kA 2 /u that is consistent with ( f4.3| ). This motion is called the Stokes drift. 

We also see why there is no net Newtonian momentum associated with phonons in a 
crystal. The atomic displacements in a harmonic crystal are given by 

r] n = Acos{k{na) — ut), (4.8) 



so the crystal equivalent of ( |4.4j) is 



^=Acos(k(na)-ut). (4.9) 

Thus 1] does not appear on the right-hand side of this equation. It is a linear equation and 
gives rise to no net particle drift. 

So, a sound wave does have real momentum? — But wait! The momentum density we 
have computed is second order in the amplitude A. The wave equation we have used to 
compute it is accurate only to first order in A. We may expand the velocity field as 

v = v (1) +v (2) + ••-, (4.10) 

where the second-order correction V( 2 ) arises because the equations of fluid motion are non- 
linear. This correction will possess both oscillating and steady components. The oscillatory 
components arise because a strictly harmonic wave with frequency ujq will gradually develop 
higher frequency components due to the progressive distortion of the wave as it propagates. 
(A plane wave eventually degenerates into a sequence of shocks.) These distortions are 
usually not significant in considerations of energy and momentum balance. The steady 
terms, however, represent 0(A 2 ) alterations to the mean flow caused by the sound waves, 
and these may possess energy and momentum comparable to that of the sound field. 
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Fig 7. Momentum Rux without mass-How. 

For example, we may drive a transducer so as to produce a beam of sound which totally 
fills a closed cylinder of fluid. At the far end of the cylinder a second transducer with 
suitably adjusted amplitude and phase absorbs the beam without reflection. Since the 
container is not going anywhere, it is clear that the average velocity of the center of mass 
of the fluid must be zero, despite the presence of the sound wave. An exact solution of 
the non-linear equation of motion for the fluid must provide a steady component in V( 2 ) 
and this counterflow completely cancels the (pivm) term. Indeed in Lagrangian coordinates 
the fluid particles simply oscillate back and forth with no net drift. The wave momentum 
and its cancelling counterflow are simply artifacts of our Eulerian description. This is one 
reason why professional fluid mechanics dislike the notion of momentum being associated 
with sound waves. 

The outlook for the two-fluid model is not entirely bleak, however. The V( 2 ) corrections 
do not always exactly cancel the momentum. Any non-zero value for V • piVti) — such as 
occurs at the transducers at the ends of the cylinder — will act as a source or sink for a 
V(2) counterflow, but its exact form depends on the shape of the container and other effects 
extrinsic to the sound field. For transducers immersed in an infinite volume of fluid, for 
example, the counterflow will take the form of a source-sink dipole field, and, although the 
total momentum of the fluid will remain zero, there will be a non-zero momentum density. 
Because it is not directly associated with the sound field, in the language of the two-fluid 
model the induced V(2) counterflow is not counted as belonging to the normal component of 
the fluid, i.e. to the phonons, but is lumped into the background superflow. It is determined 
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by enforcing mass conservation, 

V • p s v, + V • p„v n = 0. (4.11) 

The difference of opinion between the physics and fluid mechanics communities over whether 
phonons have real momentum reduces, therefore, to one of different accounting conventions. 



V. THE UNRUH WAVE EQUATION 

Now we will examine the energy and momentum in a moving fluid. The flow of an 



irrotational fluid is derivable from the action 12 



S = J d'x {p0 + l - P {V<P) 2 + u(p)} . (5.1) 

Here p is the mass density, <fi the velocity potential, and the overdot denotes differentiation 
with respect to time. The function u may be identified with the internal energy density. 
Varying with respect to <p yields the continuity equation 

p + V-(pv) = 0, (5.2) 

where v = V0. Varying p gives a form of Bernoulli's equation 

0+I v 2 + p(p) = O, (5.3) 

where p(p) = du/dp. In most applications p would be identified with the specific enthalpy. 
For a superfluid condensate the entropy density, s, is identically zero and p is the local 
chemical potential. 

The gradient of the Bernoulli equation is Euler's equation of motion for the fluid. Com- 
bining this with the continuity equation yields a momentum conservation law 

dt(pvi) + djipvjVi) + pdip = 0. (5.4) 

We simplify ( |5.4| ) by introducing the pressure, P, which is related to p by -P(p) = J pdp. 
Then we can write 
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dt(pvi) + djllji = 0, (5.5) 

where Uij is given by 

Ilij = pViVj + SijP. (5.6) 

This is the usual form of the momentum flux tensor in fluid mechanics. 

The relations p = du/dp and p = dP/dp show that P and u are related by a Legendre 
transformation: P = pp — u(p). From this and the Bernoulli equation we see that the 
pressure is equal to minus the action density: 

-P = p0+ip(V0) 2 + M ( P ). (5.7) 

Consequently, we can write 

lUj = pd^djcf) - 8ij |p0 + ip(V0) 2 + u(p)} . (5.8) 



VI. THE ACOUSTIC METRIC 

To obtain Unruh's wave equation we set 

= 00 + 01 

P = Po + Pi- (6-1) 

Here O and p define the mean flow. We assume that they obey the equations of motion. 
The quantities 0i and pi represent small amplitude perturbations. Expanding S to quadratic 
order in these perturbations gives 

S = S + Jd 4 x j Pl 0! + X - p\ + ip o (V0!) 2 + p lV • V0! j . (6.2) 

Here v = V( ) = V0 O . The speed of sound, c, is defined by 



c 2 dp 
po dp 



(6.3) 



P0 
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or more familiarly 



dP 
dp 



(6.4) 



The terms linear in the perturbations vanish because of our assumption that the zeroth-order 
variables obey the equation of motion. 



The equation of motion for p\ derived from (|6.2|) is 



Pi = -^{0i + v-V<M- (6-5) 

c z 

Since p\ occurs quadratically, we may use (|6.5| ) to eliminate it and obtain an effective action 
for the potential <p\ only 



The resultant equation of motion for <f)\ is |i3| , |14 
( d 



This can be written as[| 



(6.6) 



(6.7) 



where 



--daV^gg^d^x = 0, 



-g 



c » v, vv — c 1 



(6.8) 



(6.9) 



We find that \/—g = Jdetg^ = Pq/c, and the covariant components of the metric are 



c 



c 2 — v 2 , V T 



V, -1 

The associated space-time interval can be written 



(6.10) 



X I use the convention that Greek letters run over four space-time indices 0, 1, 2, 3 with 
Roman indices refer to the three space components. 
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= t, while 



ds 2 = ^ {c 2 dt 2 - d^dx 1 - v i dt)(dx j - v J dt)) , (6.11) 



and metrics of this form, although without the overall conformal factor po/c, appear in the 
Arnowitt-Deser-Misner (ADM) formalism of general relativity [[OJ. There, c and — v l are 
referred to as the lapse function and shift vector respectively. They serve to glue successive 
three-dimensional time slices together to form a four- dimensional space-time [|16[]. In our 
present case, provided p /c can be regarded as a constant, each 3-space is ordinary flat R 3 
equipped with the rectangular Cartesian metric glj Pace ^ = 5y — but the resultant space-time 
is in general curved, the curvature depending on the degree of inhomogeneity of the mean 
flow V. 

In the geometric acoustics limit, sound will travel along the null geodesies defined by g^ v . 
Even in the presence of spatially varying p we would expect the ray paths to depend only 
on the local values of c and v, so it is perhaps a bit surprising to see the density entering 
the expression for the Unruh metric. An overall conformal factor, however, does not affect 
null geodesies, and thus variations in po do not influence the ray tracing. 



VII. SECOND-ORDER QUANTITIES 

We are going to derive various energy and momentum conservation laws from our wave 
equation. Before we do, let us consider what sort of quantities we would want them to 
contain. 

It is reasonable to define the momentum density and the momentum flux tensor associ- 
ated with the sound field as the second order averages 

j (phonon) = (piv (1 )) + v(p 2 ), (7.1) 

and 

n (phonon) = po ( V{1) . V{1) .} _|_ Vi(piVcy)j) + Vj(piV {1)i ) + S ij (P 2 ) + V i V j (p 2 ). (7.2) 
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In these expressions I have taken no account of any steady part of V( 2 ) . This is not a quantity 
intrinsic to the sound field and has to be found by methods outside the purely acoustic. The 
other second-order quantities P 2 and p 2 can be computed in terms of first-order amplitudes. 
For P 2 we combine 

and Bernoulli's equation in the form 

A /U = -0 1 -i(V0 1 ) 2 -v-V0 1 , (7.4) 

with 

— = p fZ = ^P = l. (75) 
dp ' dp 2 dp c 2 

Expanding out and grouping terms of appropriate orders gives 

Pi = -Po(0i + v- V0O = c 2 P i, (7.6) 
which we already knew, and 

P2 = -Po^W) 2 + + y . V0O 2 . (7.7) 



We see that P 2 = y/—gL where L is the Lagrangian density for our sound wave equation. 
For a plane wave (P 2 ) = 

The second order change in the density, p 2 , may be found similarly. It is 

1 l 9 /c?lnc\ . 

p 2 = p 2 __ p . 7. 8 

& po \d\npj s 

For a plane wave the time average (P 2 ) — 0, but because p 2 contains p 2 , the time average of 
this quantity is non-zero. The resulting change of volume of the fluid, or, if the volume is 
held fixed, the resulting pressure change, is the origin of the isotropic terms in the radiation 
stress tensor discussed earlier. 



19 



VIII. CONSERVATION LAWS 



Because the linear wave equation does not have access to information about counterflows 
or second-order density changes, we will not be able to derive the real energy and momentum 
fluxes from its solution. We can still derive from the wave equation what look superficially 
like conservation laws for these quantities, however, and these laws give us insight into the 
behaviour of the solutions. The conserved quantities are of course the pseudoenergy and 
pseudomomentum. 

We begin by defining a (pseudo)-energy-momentum tensor 

= -L|^. (8.1) 

Let us recall how such a tensor comes to be associated with conservation laws. Suppose 
that we have an action S(ip,g^ v ) which is a functional of some field variables ip(x^) and the 
metric g^. If we reparameterize spacetime so that the point that had coordinates x^ is now 
denoted by x M + e M , then we have </?—></? + 5 if and — > g^ u + 5g^ u , where 

Sip = e^d^ip 

5 9tIU = D^Ey + £U„. (8.2) 

Here is the covariant derivative containing the Riemann connection compatible with the 
metric. The variation in the metric comes from the computation 

ds 2 = g^ u {x)dx^dx u 

-> 9^(x a + e a )d{x» + e»)d{x v + e v ) 

= (e a d a g flI/ + g av d^e a + g fia d l/ e a )dx t "dx u 

= (D^ v + D v e li )dxi*dx u . (8.3) 

The assembly of the terms into covariant derivatives in the last line is most easily established 
by using geodesic coordinates and the fact that 5g au is a tensor. The combination D^e u + 
D v e^ is the Lie derivative, C t g^ u , of the metric with respect to the vector field e M . 
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Since a mere re-coordinatization does not change the numerical value of the action, we 
must have 



= 5S = J tfxy/^jjUD^ + D V E») 



1 5S 



1 5S' 



(8.4) 



-g Sg^u V^g Sep 

Now the equations of motion state that S is unchanged by any variation in (p, including, a 
fortiori the change Sip = e^d^ip. Thus 

2 5S /• , , „ / 2 SS 



J d' k x y /^ge u D^ (- 



■5) 



where, in the last equality, we have integrated by parts by using the derivation property of 
the covariant derivative and the expression 

1 



-g 



-A, 



(8.6) 



for the divergence of a vector. Since (|8.5|) is true for arbitrary e fl (x), we deduce that 

2 55 N 



D, 



D^T^ = 0. 



(8.7) 



\V-gSg fa , / 

Although (|8.7|) has the appearance of a conservation law, and has useful applications 
in itself, we have not yet exploited any symmetries of the system — and it is symmetries 
that lead to conserved quantities. To derive a genuine local conservation law we need 
to assume that the metric admits a Killing vector, i] fl . This means that the particular 
reparameterization = ^ is actually an isometry of the manifold and so leaves the metric 
invariant 



AjflV = D^rj u + D u rj^ = 0. 



Combining flS.8| ) with ( |8.7| ) and using the symmetry of T^ u we find that 



(8.8) 



(8.9) 



Using (|8.6p, this can be written 



(8.10) 
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Thus the 4- vector density Q M = \f—g T^ v r] v is conventionally conserved, and 

Q = J d 3 xQ° (8.11) 

is independent of the time slice on which it is evaluated. 

Since ( |6.6| ) can be written as the usual action for a scalar field 

S 2 = J d^^g^d^xd^x, (8.12) 



we have 



= arfaerfa - <r (^g^dvfadpfa ) . (8.13) 



The derivatives with raised indices in ( |8.13|) are defined by 

1 



and 



d°<j>i = g^fa = —(0! + v • V0i), (8.14) 
Poc 



9>! = g^a^fa = — Mi + v • V0i) - c 2 d4i) ■ (8.15) 



Thus 



C 2 /£ r 



Po Vc 2 

4*- (8.16) 
Po 



The last two equalities serve as a definition of £ r and p 2 . The quantity £ r is often described 
as the acoustic energy density relative to the frame moving with the local fluid velocity [|T7 . 
It is, of course, more correctly a pseudo-energy density. 

We can express the other components of ( |S.13| ) in terms of physical quantities. We find 
that 



T*° = T°* = |^(P lV(1)i + Vi £ p )) 
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C 2 . . 

= -3 (PiU(i)i + ^P 2 ) . (8.17) 
Po 

The first line in this expression shows that, up to an overall factor, T t0 is an energy flux - 
the first term being the rate of working by a fluid element on its neighbour, and the second 
the advected energy. The second line is written so as to suggest the usual relativistic iden- 
tification of (energy- flux) /c 2 with the density of momentum. This interpretation, however, 
requires that p 2 be the second-order correction to the density, which, sadly, it is not. 
Similarly 

T 3 = — (poV(i)iV(i)j + ViPiV(i)j + Vjpiv^i + 5 ij P 2 + v i v j p 2 ) ■ (8.18) 
Po 

We again see that if we were only able to identify p 2 with p 2 then T iJ has precisely the form 
we expect for the second-order momentum flux tensor. Although it comes close, the inability 
of the pseudomomentum flux to exactly capture the true-momentum flux is inevitable as we 
know that computing the true stresses in the medium requires more information about the 
equation of state than is available to the linearized wave equation. 

We can also write the mixed co- and contra-variant components of the energy momentum 
tensor T^ v = T^ x g\ v in terms of physical quantities. This mixed tensor turns out to be more 
useful than the doubly contravariant tensor. Because we no longer enforce a symmetry 
between the indices p and u, the quantity S r is no longer required to perform double duty 
as both an energy and a density. We find 



^ T° = (E r + piV(i) • v 



cz fjT l = f ~ + v ' v (i) j (po v (i)i + PiV(p)%), (8.19) 



and 



V-gTi = -piv (1)i 

T) = - (poV(i)iV {1)j + Vipiv {l)j + <%P 2 ) • (8.20) 

We see that p 2 does not appear here, and all these terms may be identified with physical 
quantities which are reliably computed from solutions of the linearized wave equation. 
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Now we turn to the local conservation laws. In what follows I will consider only a steady 
background flow, and further one for which p , c, and hence \J—g = p\jc can be treated 
as constant. To increase the readability of some expressions I will also choose units so that 
Po and c become unity and no longer appear as overall factors in the metric or the four- 
dimensional energy-momentum tensors. I will, however, reintroduce them when they are 
required for dimensional correctness in expressions such as po v (i) or S r /c 2 . 

From the acoustic metric we find 

^o = ^(v-V)|t;| 2 

C = ~\ d i M 2 + \ V j{ d i V j - d 3 V i) 

no = \^(v-V)\v\ 2 -^d l \v\ 2 

r° = \(d iVj + d jVi ) 

r} = --Vidj \v\ 2 + -(djV k - d k Vj)(v k Vi - c 2 5 ik ) 

r) k = lv t (d jVk + d k v 3 ). (8.21) 



We now evaluate 



= d^T" + T^T^. (8.22) 



After a little algebra we find 



1 



Note the non-appearance of p\ and p 2 in the final expression — even though both quantities 
appear in T^ v . 

The conservation law therefore becomes 

d t S r + di{P lV{l)l + ViSr) + ^ij{diVj + djVi) = 0, (8.24) 



where 



^ij = PoV(i)iV(i)j + ^ijP 2 - (8.25) 
24 



This is an example of the general form of energy law derived by Longuet-Higgins and Stuart, 
originally in the context of ocean waves [E3|. The relative energy density, S r = T 00 , is not 



conserved. Instead, an observer moving with the fluid sees the waves acquiring energy from 
the mean flow at a rate given by the product of a radiation stress, S^, with the mean-flow 
rate-of-strain. This equation is sometimes cited [T9| as an explanation for the monstrous 



ship-destroying waves that may be encountered off the eastern coast of South Africa. Here 
long wavelength swell from distant Antarctic storms runs into the swift southbound Agulhas 
current and is greatly amplified by the opposing flow. 

We now examine the energy conservation law coming from the zeroth component of the 
mixed energy-momentum tensor. After a little work we find that the connection contribution 
vanishes identically and the energy conservation law is therefore 



d t (S r + piV(i) • v) + di + v ■ v (1 ))(p f(i)i + PiV(o)i)J = 0. 



(8.26) 



We see that the combination S r + piVm • v does correspond to a conserved energy — as we 
should have anticipated since a steady flow provides us with a Killing vector e = d t . This 
conservation law was originally derived by Blokhintsev for slowly varying flows, and 



more generally by Cantrell and Hart |2T[] in their study of the acoustic stability of rocket 
engines. 

Finally the covariant conservation equation D^T^ = reads 

dtp\V(i)j + di (p v {1)i v (1)j + ViPiv^j + <%P 2 ) + piV(i)idjVi = 0. (8.27) 

Here connection terms have provided a source term for the momentum density. Thus, in an 
inhomogeneous flow field, momentum is exchanged between the waves and the mean flow. 

IX. PHONONS AND CONSERVATION OF WAVE ACTION 

The conservation laws we have derived in the previous section may all be interpreted 
in terms of the semiclassical motion of phonons. As noted by Mclntyre M, the existence 
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of such an interpretation is a sure sign that the conservation laws in question are those of 
pseudomomentum and pseudoenergy. 

To make contact with the semiclassical picture we observe that when the mean flow 
varies slowly on the scale of a wavelength, the sound field can locally be approximated by a 
plane wave 

4>(x,t) = Acos(k ■ x - ujt). (9.1) 

The frequency u and the wave- vector k are related by the Doppler-shifted dispersion relation 
u) — u r + k • v, where ou r — c\k\, is the frequency relative to a frame moving with the fluid. 
A packet of such waves is governed by Hamilton's ray equations 

dx 1 duj dk l duj 



dt dk 1 ' dt dx 1 
In other words the packet moves at the group velocity 



(9.2) 



and the change in k is given by 



V g = x = + v. (9.3) 



dkj , dvi , 
—I = -k iT -^. 9.4 
dt dxo v ; 



In this equation the time derivative is taken along the ray: 

This evolution can also be derived from the equation for null geodesies of the acoustic metric 
§• 

The evolution of the amplitude A is linked with that of the relative energy density, £ r , 
through 

(Sr) = \^P^- (9-6) 

For a homogeneous stationary fluid we would expect our macroscopic plane wave to cor- 
respond to a quantum coherent state whose energy is, in terms of the (quantum) average 
phonon density N, 
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E tot = (Volume) (£ r ) = (Volume) Nhu r . (9.7) 

Since it is a density of "particles", N should remain the same when viewed from any frame. 
Consequently, the relation 

Nh=^- (9.8) 

should hold true generally. In classical fluid mechanics the quantity (£ r )/uj r is called the 
wave action |22| , p!7| . [23|| . 



The time averages of other components of the energy-momentum tensor may be also 
expressed in terms of N. For the mixed tensor we find 

(V=£ T° > = (£ r + v • Pl v (1) ) = Nhu 

(v^r ) = ((^ + v • v ( i))0boV(i)< + PiVi)) = Nfuv(V g )i 
Po 

(-V=^T°) = ( Pl v {1)i ) = Nhh 

(-V^gT)) = (p v {1 )iV (1)j + Vipiv {1)j + 5 ij P 2 ) = Nhkj(V g )i. (9.9) 

In the last equality we have used that (P2) = for a plane progressive wave. 

If we insert these expressions for the time averages into the Blokhintsev energy conser- 



vation law (|8.26|), we find that 



dt 

We can write this as 



° yhu} \- V • (NhuV g ) = 0. (9.10) 



[ ^ + V 9 ■ V«) + fuo ( ^ + V • (NV g ) ) = 0. (9.11) 



The first term is proportional to duj/dt taken along the rays and vanishes for a steady mean 
flow as a consequence of the Hamiltonian nature of the ray tracing equations. The second 
term must therefore also vanish. This vanishing represents the conservation of phonons, or, 
in classical language, the conservation of wave-action. Conservation of wave action is an 
analogue of the adiabatic invariance of E/uj in the time dependent harmonic oscillator. 
In a similar manner, the time average of ( |S.271 ) may be written 
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+ V ■ (NkjV g ) + Nki 



dvi 



) 




) 



(9.12) 



We see that the momentum law becomes equivalent to phonon-number conservation com- 
bined with the ray tracing equation (p.4|) . 



When dealing with waves in a medium it is essential to distinguish between the conceptu- 
ally distinct quantities of momentum and pseudomomentum. Under suitable circumstances 
either may be used for computing forces — but there is no general rule for determining those 
circumstances. In the case of non-dispersive sound waves, momentum and pseudomomentum 
often travel together, and are therefore easily confused — but the conservation laws derived 
by manipulating the wave equation are those of pseudomomentum and pseudoenergy. The 
GR analogy provided by the acoustic metric provides a convenient and structured route to 
deriving these laws when the background fluid is moving. 



This work was supported by NSF grant number DMR-98-17941. I would like to thank 
the staff and members of NORDITA in Copenhagen, Denmark, for their hospitality while 
this paper was being written. 



X. SUMMARY 



XI. ACKNOWLEDGEMENTS 



28 



REFERENCES 
[1] Sir Rudolf Peierls, Surprises in Theoretical Physics, (Princeton, 1979). 
[2] E. I. Blount, Bell Telephone Laboratories Technical Report 38139-9 (1971). 
[3] J. P. Gordon, Phys. Rev. A 8 (1973) 14. 
[4] D. F. Nelson Phys. Rev. A 44 (1991) 3985. 
[5] M. Stone, Phys. Rev. B 61 (2000) 11780. 

[6] G. E. Volovik, Pis'Ma Zhu. Eksp. Teor. Fiz. 67 (1998) 841 [JETP Letters 67 (1998) 
pp.881]. ( |cond-mat/9804308|) 

[7] M. E. Mclntyre, J. Fluid. Mech. 106 (1981) 331. 

[8] M. Stone, Phys. Rev. E 62 (2000) 1341. 

[9] E. J. Post, J. Acoustical Soc. Amer. 25 (1953) 55. 
[10] R. T. Beyer, J. Acoustical Soc. Amer. 63 (1978) 1025. 
[11] L. Brillouin, Annales de Physique, 4 (1925) 528 (in French). 
[12] A. M. J. Schakel, Mod. Phys. Lett. BIO (1996) 999. (|cond-mat/9607164|) 
[13] W. Unruh, Phys. Rev. Lett. 46 (1981) 1351. 
[14] W. Unruh, Phys. Rev. D 51 (1995) 282 

[15] R. Arnowit, S. Deser, C. W. Misner in Gravitation: an Introduction to Current Re- 
search ed. L. Witten (Wiley N.Y. 1962), pp227-265; C. Misner, K. Thorne, J. Wheeler, 
Gravitation (W. H. Freeman, San Francisco 1973). 

[16] For a picture see: C. Misner, K. Thorne, J. Wheeler, op. cit. p504. 

[17] Sir James Lighthill, Waves in Fluids (Cambridge University Press 1978) 

[18] M. S. Longuet-Higgins, R. W. Stuart, J. Fluid. Mech. 10 (1961) 529. 

29 



[19] M. S. Longuet-Higgins, R. W. Stuart, Deep Sea Res. 11 (1964) 529. 

[20] D. I. Blokhintsev, Acoustics of a Non-homogeneous Moving Medium. (Gostekhizdat, 
1945). [English Translation: N. A. C. A. Technical Memorandum no. 1399, (1956)] 

[21] R. H. Cantrell, R. W. Hart, J. Acoustical Soc. Amer. 36 (1964) 697. 

[22] C. J. R. Garrett, Proc. Roy. Soc. A299 (1967) 26. 

[23] D G. Andrews, M. E. Mclntyre, J. Fluid. Mech. 89 (1978) 647. 



30 



